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^ ■ ABSTRACT 

o : 

Context. Our understanding of the star formation process has traditionally been confined to certain mass or luminosity boundaries because most 
(-H Studies focus only on low-, intermediate- or high-mass star-forming regions. Therefore, the processes that regulate the formation of these different 
^ objects have not been effectively linked. As part of the "Water In Star-forming regions with HerscheF (WISH) key program, water and other 
h— 5 important molecules, such as CO and OH, have been observed in 5 1 embedded young stellar objects (YSOs). The studied sample covers a range 
of luminosities from < 1 to > 10' Lq. 

Aims. We analyse the CO line emission towards a large sample of embedded protostars in terms of both line intensities and profiles. This analysis 
covers a wide luminosity range in order to achieve a better understanding of star formation without imposing luminosity boundaries. In particular, 
this paper aims to constrain the dynamics of the environment in which YSOs form. 
' Methods. Herschel-mF\ spectra of the '^CO 7= 10-9, "CO 7= 10-9 and C'*0 7=5-4, 7=9-8 and 7=10-9 lines ai'e analysed for a sample of 5 1 
embedded protostars. In addition, JCMT spectra of '^CO 7=3-2 and C'*0 7=3-2 extend this analysis to cooler gas components. We focus on 
\^ characterising the shape and intensity of the CO emission line profiles by fitting the lines with one or two Gaussian profiles. We compare the 
values and results of these fits across the entire luminosity range covered by WISH observations. The effects of different physical parameters as a 
function of luminosity and the dynamics of the envelope-outflow system are investigated. 

Resuhs. All observed CO and isotopologue spectra show a strong linear correlation between the logarithms of the line and bolometric luminosities 
^ across six orders of magnitude on both axes. This suggests that the high-7 CO lines primarily trace the amount of dense gas associated with 
$^ \ YSOs and that this relation can be extended to larger (extragalactic) scales. The majority of the detected '"CO line profiles can be decomposed 
, into a broad and a narrow Gaussian component, while the C'^O spectra are mainly fitted with a single Gaussian. For low- and intermediate-mass 
■ protostars, the width of the C'^O 7=9-8 line is roughly twice that of the C'^O 7=3-2 line, suggesting increased turbulence/infall in the warmer 
I I inner envelope. For high-mass protostars, the line widths are comparable for lower- and higher-7 lines. A broadening of the line profile is also 
observed from pre-stellar cores to embedded protostars, which is due mostly to non-theiTnal motions (turbulence/infall). The widths of the broad 
'-CO 7=3-2 and 7= 10-9 velocity components correlate with those of the narrow C '**0 7=9-8 profiles, suggesting that the entrained outflowing gas 
^ and envelope motions are related independent of the mass of the protostar. These results indicate that physical processes in protostellar envelopes 
OO have similar characteristics across the studied luminosity range. 

in 
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1. Introduction stars (e.g.. fShepherd & Church welll 1 1996t iBeuther & Shepherd 

m , • . • , , , J , • • • , 2005). The reason is related to the short life-time (IMottram et al 



The evolution of a protostar is closely related to the initial mass 201 1) and the large distances (few kpc) associated with massive 
of the molecular core from which it forms and to the spe- ysOs. This means that outflows from massive stars are less well 
cific physical an d chemica properties of the original molec - ^^^^j^^^ ^^^^ ^^^^^ low-mass counteipaits. The agent that drives 



• i-H , , , / ITT; — 1 1 ' f I — ' ' — I I I icsuivcu Liiaii Liicii luw-mass cuuiiLcinaiLS. 1 11c agciiL Liiai uiivcs 

V ■ i ^i^';,^^"'!^^^-^.; I^^^^lf 1^ | vanDishoeck& Blake| | l998t molec ular outflow, either jets or win ds from the disk and/or 

5—1 < . 7, — , . . . . . . . . stars (e.g.. lChurchwelllll999HArce et al.ll2007i) . might be differ- 

, tion, young ste lar objects (YSOs) are embedded in large cold depending on the mass of the star-forming region. Therefore, 

■ ■ ■ and dusty envelopes which will be accreted or removed by the ^j^^ interaction of the outflow with the suiTounding material and 

forming star. Depending on the niass of the star-forming region, ^ i^jj ^^e resulting chemisny may differ across the mass 

the parameters and mechanisms that rule several processes of the j-^jjee 

star formation, such as driving agent of the molecular outflow ^, . , ,.rr , 

and accretion rates, will vaiy. The accretion rates are also different depending on the mass 

Molecular outflows are crucial for removing angular mo- *^ ^°™7"^, '^f ; J^P^'^f ^ o!^""?/"' ^"^""^^^^ ^^^i", f°™|r 

mentum and mass from the protostellar system (see review by ^1°" '^K ^ ^ ' ^^^"^^^ 

ra^999l) . Thev have been extensi velv studied for low -mass ^^ereas higher values are necessary in order to overcome ra- 

YSOs (e.g., ICabrit &Bertoud [T9921 iBachiller & Tafalli \l99^ diatio n pressure and fomi m assive stars within a free-fall time 

where they are better characterised than for massive proto- e-gj. Unma & Ada^ M- These values range from 10"^ to 

10 Moyr ' for sources with > 10 Lq (e.g.. IBeuther et al.l 

• Herschel is an ESA space observatory with science instruments pro- 2002). In addition, for the low-mass sources, the accretion 

vided by European-led Principal Investigator consortia and with impor- episode finishes before the protostar reaches the main sequence, 

tant participation from NASA. while massive YSOs still accrete circumstellar material after 
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reaching the h ydrogen burning phase (iPalla & Stahlerl Il993i 
ICesaronill2005h . 

Another difference is that the ionising radiation created by 
main-sequence OB stars is much more powerful than that gen- 
erated by a single low-mass protostar. Therefore, photon domi- 
nated regions (PDR) and Hii regions are formed in areas of mas- 
sive star formation, affecting the kinema tics, temperature and 
chemistry of the surrounding material ( HoUenbach & TielensI 
[1991 . In addition, the strength of stellar winds and their inter- 
action with the envelope material is different depending on the 
stellar spectral type of the YSOs. 

Because of these differences, the study of star formation has 
traditionally been restricted to mass boundaries, focused either 
on low-mass (M < 3 Mq) or high-mass (M > 8 Mq) YSOs. 
One of the goals of the "Water In Star-forming regions with 
Herschel" (WISH) key program (van Dishoeck et al. 2011) is 
to offer a complete description of the interaction of young stars 
with their surroundings as a function of mass. For this purpose, 
and in order to constrain the physical and chemical processes 
that determine star formation, water and other key molecules 
like CO have been observed for a large sample of embedded 
YSOs (51 sources). The targeted objects cover a vast range of 
luminosities (from < 1 L© to > 10^ Lq) and different evolution- 
ary stages (more details in Section 12.111. With the Heterodyne 
Instrument for the Far-Infrared (HIFI; 'de Graauw et al. 2010) 
on board the Herschel Space Observatory (Pilbratt et al. 2010), 
high spectral resolution data of high frequency molecular lines 
have been obtained. These can be used to probe the physical con- 
ditions, ch emical composition and dynamics of protostellar sys- 
tem s (e.g..'Evans'1 999l;|j0rgensen et al.l2002tlvan der Wiel et all 
12013. s ubmitted). 

Due to its high, stable abundance and strong lines, CO is 
one of the most important and often used molecules to probe 
the different physical components of the YSO environment (en- 
velope, outflow, disk). In particular, molecular outflows are 
traced by '^CO e mission through maps in the line wings (e.g, 
ICurtis et aLOOlOl) . Its isotopologue C'**0 is generally thought to 
probe quiescent gas in the denser part of the protostellar enve- 
lope, whereas '^CO lines origi nate in the extended envelope and 
the outflow cavity walls (e.g.. ISpaans et al.lll995l : lOraves et al.l 
l2010t lYildiz et al.ll2012h . In addition, CO has a relatively low 
critical densities, due to its small permanent dipole moment 
(~0.1 Debye), and relatively low rotational energy levels, so this 
molecule is easily excited and thermalised by collisions with H2 
in a typical star formation environment. For this reason, mea- 
surements of CO excitation provide a trustworthy estimate of the 
gas kinetic temperature. Moreover, integrated intensity measure- 
ments can be used to obtain column densities of warm gas, pro- 
viding a reference to determine the abundances of other species, 
such as water and H2. 

Most CO observations from ground-based sub-millimetre 
telescopes have been limited to low-7 rotational transitions (up 
to upper transition 7u=3, i.e., upper-level energy E^jk^ ~35 K), 
or mid-y transitions { Ja-(), with a E^/k^ of ~100 K). Thanks to 
HIFI, spectrally resolved data for high-7 CO transitions (J^ up 
to 16, Eti/kB ~600 K) are observable for the first time, so warm 
gas di rectly associated with the forming star is probed (e.g., 
Yildiz et al. 2010, 2012: iPlume et all 120121: Ivan der Wiel et all 
2013. submitted) . Therefore, a uniform probe of the YSOs over 
the entire relevant range of E^/kB (from 10-600 K) is achieved 
by combining HIFI data with complementary spectra from 
single-dish ground-based telescopes. These observations are in- 
dispensable in order to ensure a self-consistent data set for anal- 
ysis. Finally, the study of these lines in our Galaxy is crucial 



in order to compare them with the equivalent lines targeted in 
high-redshift galaxies which are often used to determine star- 
formation rates on larger scales. 

In this paper we present '^CO 7=10-9, '^CO 7=10-9, 
C^**0 7=5-4, 7=9-8 and 7=10-9 HIFI specti-a of 51 YSOs. 
Complementing these data, '^CO and C"^0 7=3-2 spectra ob- 
served with the James Clerk Maxwell Telescope (JCMT) are in- 
cluded in the analysis in order to use CO to its full diagnostic 
potential and extend the analysis to different regions of the proto- 
stellar environment with different physical conditions. Section|2] 
describes the sample, the observed CO data and the method de- 
veloped to analyse the line profiles. A description of the mor- 
phology of the spectra, an estimation of the kinetic temperatures 
and correlations regarding the line luminosities of each isotopo- 
logue transition are presented in Section|3] These results are also 
compared to other YSO parameters such as luminosity and en- 
velope mass. In Section|4]we discuss the results to constrain the 
dynamics of individual velocity components of protostellar en- 
velopes, characterise the turbulence in the envelope-outflow sys- 
tem and consider high-7 CO as a dense gas tracer Our conclu- 
sions are summarised in Section |5] 



2. Observations 

2.1. Sample 

The sample discussed in this paper is drawn from the WISH 
survey and covers a wide range of luminosities and different 
evolutionary stages. A total of 51 sources are included in this 
study, which can be classified into three groups according to 
their bolometric luminosities, Lboi- The sub-sample of low-mass 
YSOs, characterized by Lboi < 50 Lp, is composed of 15 Class 
and 1 1 Class I protostars (see lEvans et aL 1 120091 for details of 
the classification). Six intermediate-mass sources were observed 
with 70 Lq < Lboi < 2x10^ Lq. Finally, 19 high-mass YSOs 
with Lboi > 2x10^ L0 complete the sample. The bolometric lu- 
minosity of the sample members, together with their envelope 
masses (Menv), distances (d) and source velocities (ulsr) are 
summarised in Table [1] For more information about the sam- 
ple studied in WISH, see van Dishoeck et al. (2011). Focusing 
on the evolutionary stages, the sub-sample of low-mass YSOs 
ranges from Class to Class I, the intermediate-mass objects 
from Class to Class I as well, and in the case of the high-mass 
sources, from (mid-IR-quiet/mid-IR-bright) massive young stel- 
lar objects (MYSOs) to ultra-compact H n regions (UCH 11). 

2.2. HIFI observations 

The sources were observed with the Heterodyne Instrument for 
the Far Infrared (HIFI) on the Herschel Space Observatory. 
The HIFI CO and isotopologue lines studied in this paper are: 
'^CO 7=10-9, '^CO 7=10-9, C'**0 7=5-4, 7=9-8 and 7=10-9. 
The upper-level energies and frequencies of these lines together 
with the HIFI bands, main beam efficiencies (t/mb), beam sizes, 
spectral resolution and integration times are presented in Table 
|2] With the exception of the '^CO 7=10-9 line, all isotopologue 
line observations were obtained together with H2O lines. The 
'^CO 7=10-9 line was targeted for the low- and intermediate- 
mass sample but only for one high-mass object (W3-IRS5). 
The i^CO 7=10-9 and C''^0 7=9-8 Hnes were observed for 
the entire sample, while C'^O 7=5-4 only for the Class and 
intermediate-mass protostars. C'^O 7=10-9 was observed for all 
low-mass Class sources, two low-mass Class I (Elias 29 and 
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Table 1. Source parameters. 



Source 


fLSR 




d 


-^env 


References 




(km s-i) 


(Le) 


(kpc) 


(Mo) 




Low-mass: Class 












L 1448-MM 


-^5.2 


9.0 


0.235 


3.9 


1 


NGC 1333 IRAS 2 A 


+1.1 


35.7 


0.235 


5.1 


1 


NGC 1333 IRAS4A 


+1.Q 


9.1 


0.235 


5.6 


1 


NGC 1333 IRAS4B 


+1.4 


4.4 


0.235 


3.0 


1 


L1527 


+5.9 


1.9 


0.140 


0.9 


1 


CedllOIRS4 


+4.2 


0.8 


0.125 


0.2 


1 


BHR71 


-4.4 


14.8 


0.200 


2.7 


1 


IRAS 15398 


+5.1 


1.6 


0.130 


0.5 


1 


L483-MM 


+5.2 


10.2 


0.200 


4.4 


1 


Ser SMM 1 


+S.5 


30.4 


0.230 


16.1 


1 


SerSMM4 


+S.Q 


1.9 


0.230 


2.1 


1 


Ser SMM 3 


+1.6 


5.1 


0.230 


3.2 


1 


L723-MM 


+ 11.2 


3.6 


0.300 


1.3 


1 


B335 


-^8.4 


3.3 


0.250 


1.2 


1 


LI157 


+2.6 


4.7 


0.325 


1.5 


1 


Low-mass: Class I 












LI489 


+7.2 


3.8 


0.140 


0.2 


1 


LI551 IRS 5 


+6.2 


22.1 


0.140 


2.3 


1 


TMRl 


+6.3 


3.8 


0.140 


0.2 


1 


TMC lA 


+6.6 


2.7 


0.140 


0.2 


1 


TMC 1 


+5.2 


0.9 


0.140 


0.2 


1 


HH46 


+5.2 


27.9 


0.450 


4.4 


1 


IRAS 12496 


+3.1 


35.4 


0.178 


0.8 


1 


Ellas 29 


+4.3 


14.1 


0.125 


0.3 


1 


Oph IRS 63 


+2.8 


1.0 


0.125 


0.3 


1 


GSS30IRSI 


+3.5 


13.9 


0.125 


0.6 


I 


RN0 91 


+0.5 


2.6 


0.125 


0.5 


I 


Intermediate-mass 












NGC 7 129 FIRS 2 


-9.8 


430 


1.25 


50.0 


2 


L1641 S3MMS1 


5.3 


70 


0.50 


20.9 


2 


NGC 2071 


9.6 


520 


0.45 


30.0 


2 


Vela IRS 17 


3.9 


715 


0.70 


6.4 


2 


Vela IRS 19 


12.2 


776 


0.70 


3.5 


2 


AFGL490 


-13.5 


2000 


1.00 


45.0 


2 


High-mass 












IRAS05358-I-3543 


-17.6 


6.3 X lO'' 


1.8 


142 


3 


IRAS 16272-4837 


-46.2 


2.4 X 10"* 


3.4 


2170 


3 


NGC6334-I(N) 


-3.3 


1.9 X 10^ 


1.7 


3826 


3 


W43-MM1 


+98.8 


2.3 X lO'* 


5.5 


7550 


3 


DR21(0H) 


-3.1 


1.3 X 10"* 


1.5 


472 


3 


W3-IRS5 


-38.4 


1.7 X 10' 


2.0 


424 


3 


IRAS 18089- 1732 


+33.8 


1.3 X lO'* 


2.3 


172 


3 


W33A 


+37.5 


I.I X 10' 


3.8 


1220 


3 


IRAS18151-1208 


+32.8 


2.0 X lO'* 


3.0 


153 


3 


AFGL259I 


-5.5 


2.2 X 10' 


3.3 


320 


3 


G327-0.6 


-45.3 


5.0 X 10-* 


3.3 


2044 


3 


NGC6334-I 


-7.4 


2.6 X 10' 


1.7 


500 


3 


G29.96-0.02 


+97.6 


3.5 X 10' 


6.0 


768 


3 


G3I.41-hO.31 


+97.4 


2.3 X 10' 


7.9 


2968 


3 


G5.89-0.39 


+ 10.0 


5.1 X lO* 


1.3 


140 


3 


GlO.47+0.03 


+67.3 


3.7 X 10' 


5.8 


1168 


3 


G34.26+0.I5 


+58.0 


3.2 X 10' 


3.3 


1792 


3 


W51N-eI 


+59.5 


I. Ox 10' 


5.1 


4530 


3 


NGC7538-IRS1 


-56.2 


1.3 X 10' 


2.7 


433 


3 



Notes. See lvan Dishoeck et alj ( 1201 If) for the source coordinates. 

References. (1 ) Bolometric l uminosities and envelope masses obtained from iKristensen et al.l l l2012h . (2) Envelope ma sses collected in 
IWampfler et al.l (I20I3. in press!) . (3) Bolometric luminosities (obtained from observations) and envelope masses calculated in Ivan der Tak et si] 
||20I3. submittecl)" 
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GSS 30 IRS 1), one intei-mediate-mass YSO (NGC 7129) and the 
entire high-mass sub-sample. 

Single-pointing observations were performed for all targets 
in dual-beam-switch (DBS) mode, chopping to a reference posi- 
tion 3' from the target. There is no contamination due to emis- 
sion at the off position except for the '^ CO 7=10-9 spectr um 
of NGC1333IRAS2A and IRAS4A (see lYildiz etlDlioiOl for 
more details). These spectra have been corrected and presented 
in this paper without contamination. In the case of W43-MM1, 
the absorption fe atures found in the '^C O 7=10-9 spectrum are 
caused by H2O+ dWvrowski et al .1120101) . 

HIFI has two backends: the Wide Band Spectrometer (WBS) 
and the High Resolution Spectrometer (HRS). Both spectrome- 
ters simultaneously measure two pola rizations, horizontal (H ) 
and vertical (V). For more details, see 'R oelfsema et al.l (1201 2h . 
The WBS has a constant spectral resolution of 1 . 1 MHz, whereas 
the HRS has different configuration modes with four possible 
spectral resolutions: 0.125, 0.25, 0.5 and 1.0 MHz. The spec- 
tral resolution for each of the studied HIFI lines is listed in 
Table |2] The WBS data present lower noise than the HRS data 
(factor of V2) and provide a good compromise between noise 
and resolution. Therefore, the WBS data are the primary focus 
of this paper HRS observations are only used for analysing the 
C'^O 7=5^ line for the low-mass sources because their narrow 
line profiles require the higher spectral resolution provided by 
these data. 

The data reduction was performed using the standard HIFI 
pipeline in the Herschel Interactive Processing Environment 
(HIPli3) ver. 8.2 (Ott 2010), resulting in absolute calibration 
on the corrected antenna temperature scale, and velocity 
calibration with a dlsr precision of a few ms"'. The version 
of the calibration files used is 8.0, released in February 2012. 
The flux scale accuracy was estimated to be 10% for bands 
1, 4 and 5. Subsequently, the data were exported to GILDAS- 
CLAS^ for further analysis. The H and V polarizations were 
observed simultaneously and the spectra averaged to improve 
the signal-to-noise ratio. In order to avoid possible discrep- 
ancies between both signals, the two polarisations were in- 
spected for all the spectra presented in this paper with no dif- 
ferences > 20% found. Afterwards, line intensities were con- 
verted to main-beam brightness temperatures through the rela- 
tion Tmb - /?7mb (see Wilson et al. 2009 for further informa- 
tion about radio-astronomy terminology). The main beam effi- 
ciency , 77MB, for each HIFI band was taken from lRoelfsema et al.l 
( I2OI2I) and listed in Table|2] The final step of the basic reduction 
was the subtraction of a constant or linear baseline. 

2.3. JCMT ground-based observations 

Complementary data from the 15-m James Clerk Maxwell 
Telescope (JCMT) on Mauna Kea, Hawaii are also included in 
this paper, in particular for the high-mass sources for which 
'■^CO 7=10-9 data are not available. Jiggle map observations 
of i^co 7=3-2 and C'-^O 7=3-2 for a sub-sample of YSOs 
were obtained ^ith the Heterodyne Array Receiver Program 
(HARP Buc kle etalJ l2009) in August 2011 and summer 2012 
(proposal Ml 1BN07 and M12BN06). For the sources and tran- 
sitions not included in the proposal, comparable data were ob- 
tained from the JCMT public archive. Four low-mass sources 

' HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Centre, and the HIFI, PACS and SPIRE consortia. 

2 http://www.iram.fr/IRAMFR/GILDAS/ 



were observed with the 12-m Atacama Pathfinder Experiment 
Telescope, APEX, because these protostars are not visible 
from the JCMT (see Appendix |B]). Further inforn iation about 
the low-mass Y SOs and data can be found in lYildiz et al.l 
(I2OI3. submittedh . 

The HARP instrument is a 4x4 pixel receiver array, although 
during the observation period one of the receivers (H14) was 
not operational. The lines were observed in position-switching 
mode, with the off-positions carefully chosen in order to avoid 
contamination. For the most massive and crowded regions, test 
observations of the off-position were taken for this purpose. The 
spatial resolution of the JCMT at the observed frequencies is 
~14", with a main beam efficiency of 0.6fl This same value of 
T/MB was used for the data obtained from the JCMT archive be- 
cause the small variations of this parameter (< 10%) recorded 
over time are negligible compar ed to the calibration uncertain- 
ties of the JCMT ( -20% . iBuckie' etal. 2009). Some of the spec- 
tra collected from the JCMT archive were observed in a lower 
spectral resolution setting. Therefore, for these data the spec- 
tral resolution is 0.4 kms"' instead of 0.1 kms"' (indicated in 
TableEJ. 

In the first step of the reduction process, the raw ACSIS data 
downloaded from the JCMT archive were transformed from sdf 
format to fits format using the StarlinlsQ package for each and ev- 
ery pixel. Next, the data were converted to CLASS format and the 
central spectrum was extracted after convolving the map to the 
same beam size as ffie ^^CO 7=10-9 HIFI observations (20")- 
Line intensities were then converted to the main-beam bright- 
ness temperature scale and linear baselines subtracted. Since this 
manuscript focuses on analysing and comparing the central spec- 
trum of the studied YSOs, the full JCMT spectral maps will be 
presented and discussed in a forthcoming paper 

2.4. Decomposition method 

In order to quantify the parameters that fit each spectrum, the fol- 
lowing procedure was applied to all spectra. First, the data were 
re-sampled to 0.27 km s ' so that the results could be compared 
in a systematic manner. Then, the spectra were fitted with a sin- 
gle Gaussian profile using the IDL function mpfitfun, after which 
we plotted the residuals obtained from the fit to confirm whether 
the line profile hid an additional Gaussian component. For those 
sources whose profiles showed clear sub-structure, i.e., the resid- 
uals were larger than 3 sigma rms, a two Gaussian component fit 
was used instead. Examples of the decomposition procedure are 
shown in Fig.[T] The results of this process together with the rms 
and integrated intensity for all lines are presented in Tables [ATI 
to lA.5l in AppendixE] 

All HIFI lines are observed in emission and none of the 
HIFI spectra present clear infall signatures. Moreover, some CO 
lines show weak self-absorption features, which are of marginal 
significance and will not be discussed further In addition, ex- 
tremely high velocity (EHV) emission features have been iden- 
tified. The EHV components are knotty structures spaced regu- 
larly and associated with sh ocked material movin g at velocities 
of hundreds of kms ' (e.g. iBachiller et aT]ll990h . These struc- 
tures have been detected in the '^CO 7=10-9 spectra for the low- 
mass Class_0_sou^ BHR71 ( Kristensen et al] 
I20T1I lYildiz et al.l I2OI3. submitted!) ■ These EHV components 
were not included in the study of the line profiles, so the resid- 



^ http://www.jach.hawaii.edu/JCMT/spectral_line/General/status.html 
* http://starlink.jach.hawaii.edu/starlink 
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Table 2. Overview of the main properties of the observed lines. 



Mol. 


Trans. 


Eu/ks 


Frequency 


Tel./Inst.-band 


'7MB 


Beam 


Spec. Resol. 


Exposure time (min) 






(K) 


(GHz) 






size (") 


(kms"') 


LMO 


LMI 


IM 


HM 




3-2 


33.2 


345.796 


JCMT 


0.63 


14 


0.1/0.4 
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Notes. LMO: low-mass Class sources. LMI: low-mass Class I objects. IM: intermediate-mass YSOs. HM: high-mass protostars. 
<"' WBS data. HRS data. NGC1333 IRAS 2A, 4A and 4B were observed for 300 min each. 



uals were analysed after fitting each of these features with a 
Gaussian function and subtracting them from the initial profile. 

The meth od used for exam ining the data is similar to that 
introduced by iKristensen et al.l (|201 Oi) for several water lines in 
some of the WISH low-mass YSOs, applied to high-7 CO in 
lYildizet all fcolO) and extended in Kristensen et al. (2012) for 
the 557 GHz lio-loi water line profiles of the entire low-mass 
sample. The emission lines are classified as narrow or broad if 
the full width half maximum (FWHM) of the Gaussian compo- 
nent is lower or higher than 7.5 kms"', respectively. This dis- 
tinction is made because 7.5 kms"' is the maximum width ob- 
tained in the single Gaussian fit of the HIFI C'^O lines, which is 
considered as naiTow and traces the dense warm quiescent enve- 
lope material (see Section|4]for further analysis and discussion). 

The naiTow component identified in the high-/ CO isotopo- 
logue lines is always seen in emission, unlike a co mponent of 
simila r width seen in the lio-loi line of water ( Kris tensen et al.l 
l2012h . The narrow components in high-7 CO and low-7 water 
probe entirely different parts of the protostar: the former traces 
the quiescent warm envelope material, the latter traces the cold 
outer envelope and ambient cloud. The broad emission in low- 
J CO is typically narrower than in water and traces entrained 
outflow material. Only the highest-/ lines observed by HIFI 
trace the same warm shocke d gas as seen in the water lines 
dYildiz et all 1201 3. submittedh . To summarise, the components 
identified in the CO and isotopologue data cannot be directly 
compared to those observed in the H2O 1 lo-loi lines because the 
physical and chemical condi tions probed by water are different 
to tho se probed by CO (see ISantangelo et 311120121: IVasta et aTl 
I2012h . 

In the analysis of the JCMT data, the FWHM of the broad 
velocity component for the complex '^CO J-3-2 line profiles 
was disentangled by masking in each spectrum the narrower 
emission and self-absorption features. The width of the nar- 
row C'^O 7=3-2 lines was constrained by fitting these profiles 
with a single Gaussian. The results of these fits are presented in 
TablelEUin AppendixE] 



3. Results 

One of the aims of this paper is to characterise how the observed 
emission lines compare as a function of source luminosity. In 
order to simplify the comparison across the studied mass range, 
the main properties and parameters of the HIFI and JCMT lines, 
such as line morphology, total intensity and kinetic temperature, 
are presented in this section. A more detailed description of the 



line profiles is reserved for Appendices lAland |B](Figs . lA. 1 I to IA.5I 
show the HIFI spectra and Figs.|B3]and|B3]the JCMT data). 

Further study and analysis of each sub-sample will be pre- 
sented in several forthcoming papers. The CO lines for the 
low-mass sources and their excitation will be discussed by 
Yildiz et al. (2013, submitted). A review of the intermediate- 
mass sources focused o n the water lines will be performed by 
iM'^Coev et al.l din prep.l) . In the case of high- mass YSOs, low-7 
H2O line profiles will be studied in detail by Ivan der Tak et al.l 
(12013. submitted!) ■ 

In this manuscript a summary with the main characteristics 
of the studied emission line profiles is presented in Section [TTI 
Section U!2l describes the calculation of the line luminosity, Leo, 
for each observed isotopologue together with its correlation with 
Lboi- Finally, in Section |331 an estimation of the kinetic temper- 
ature is obtained for two sources, an intermediate-mass and a 
high-mass YSO, and compared with values obtained for low- 
mass sources. 

3. 1 . Characterisation of tiie line profiles 

Figures |2] and [3] show characteristic profiles of each transition 
and YSO sub-type, so the line profiles can be compared across 
the luminosity range. '^CO 7=10-9 spectra present more intense 
emission lines than the other observed isotopologues and more 
complex line profiles. Two velocity components are identified 
and most of the '^CO 7=10-9 spectra can be well fitted by two 
Gaussian profiles (Fig. [T^). Weak self-absoiption features are 
also observed in some sources, such as SerSMMl. '^CO 7=10- 
9 profiles are weaker and narrower than '^CO 7=10-9 spectra. 
Some of the detected lines, especially for the high-mass sam- 
ple, are fitted using two Gaussian components (Fig. [T^). In the 
case of the C"*0 7=5-4 spectra, a weak broad velocity compo- 
nent is identified in 6 sources (indicated in Table IA.3b . due to 
the long exposure time and the high S /N reached for this transi- 
tion. The width of this broad component is narrower by a factor 
of 2-3 than that detected for the '^CO and '^CO 7=10-9 lines. 
Similarly, two velocity components have been identified for 
file C'^O 7=9-8 line in three high-mass sources: GlO.47+0.03, 
W51N-el and G5.89-0.39 (see Fig.lA3]). These massive objects 
present the widest broad velocity components for both '^CO and 
C'^^O 7=10-9 spectra. The width of the broad C"^0 7=9-8 com- 
ponent is slightly smaller than that identified in the '^CO 7= 10-9 
emission for each of these YSOs. 

Two velocity components were previously identified in ap- 
proximately half of the 20 deeply en ibedded young stars in the 
Taurus molecular cloud studied by iFuUer & Laddl (l2002l) us- 
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Fig. 1. Gaussian decomposition for the CO and isotopologues line profiles: a) two Gaussian fit for the line profiles with two different 
velocity components identified, such as the '^CO 7=10-9 spectra for the low-mass YSOs SerSMM3 (left), the intermediate-mass 
Vela 19 (centre) and the '^CO 7=10-9 spectrum of the high-mass W33A (right), b) Single Gaussian fit of sources characterised 
by one component profile, such as the C'^O 7=9-8 spectra of W33A. The red lines show the Gaussian fits and the green lines the 
baseline. 
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Fig. 2. i^co 7=10-9 (left) and '^CO 7=10-9 (right) spectra 
for a low-mass Class protostar (top, SerSMMl), low-mass 
Class I source (Elias 29), intermediate-mass object (NGC2071) 
and high-mass YSO (bottom, W3-IRS5). The green line indi- 
cates the baseline level and the red dashed line the kms ' 
value. All spectra have been re-binned to 0.27 km s"' and shifted 
with respect to their relative local standard of rest velocity. 



ing lower-7 C'^O observations. They found typical FWHM line 
widths of ~0.6 and ~2.0 kms ' for the narrow and broad com- 
ponent respectively. These values are significantly narrower than 
the widths obtained from the HlFl data, so our interpretation and 
analysis of these com ponents is different to that presented by 
iFuUer & Laddl (|2002|) . On the other hand, the bulk of C'^^O line 
profiles (especially 7=9-8 and 7=10-9 transitions) are generally 
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Fig. 3. Same as Fig.|2]but for the C'**0 spectra from the observed 
transitions: 7=5^ (left), 7=9-8 (centre) and 7=10-9 (right). For 
details about these objects see AppendixlAl 



well fitted by a single Gaussian (for an example, see Fig. [TJ)). 
Therefore, in our analysis only the narrow C"^0 component is 
considered for the three sources with two velocity component 
profile. Regarding the line intensity, the spectra of the observed 
high-mass YSOs have higher main beam temperatures than the 
spectra of the intermediate-mass objects, which in turn show 
stronger lines than the low-mass sources. 

Another result obtained when we extend this characterisa- 
tion to the JCMT data is the complexity of the '^CO 7 =3-2 
profiles compared to the '^CO 7=10-9 spectra (see Figs. IA.11 
and IB. 41 for comparison across the studied sample). The HIFI 
data probe warmer gas from inner regions of the molecular core 
and present simpler emission line profiles (with no deep ab- 
sorptions and foreground emission features) than the lower-7 
spectra. However, similar to the '^CO 7=10-9 line profiles, the 
'^CO 7=3-2 spectra can be decomposed into different velocity 
components. A broad velocity component is identified in 39 out 
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of 47 sources, ranging from ~7.4 to 53.5 kms"' in width. For 
C'^O, the shape of the 7=3-2 lines are very similar to those of 
the 7=9-8 lines (see Figs. lB.ll to lB.3l for examples). 

The FWHM of the ^^CO 7=3-2 broad component for most 
of the high-mass sources is approximately double the width 
obtained for the '^CO 7=10-9 broad component (values in 
Tables IASI and lBTTT i. In the case of the one source for which a 
'^CO 7=10-9 observation was performed as part of WISH (W3- 
IRS5), the width of the broad component is similar to that calcu- 
lated for the 7=3-2 spectrum (factor of 1.2±0.1) and is twice the 
width of the '^CO 7=10-9 emission line (see Fig.[2ll. Sim ilar ra- 
tios were found bv lvan der Wiel et al. (2013. submitted) for the 
high-mass source AFGL2591 as part of the CHESS ("Chemical 
HErschel Survey of Star-forming regions") key program ob- 
servations. For the intermediate-mass object NGC2071, the 
'^CO 7=10-9 broad component is 1.7±0.1 larger than the width 
of the '^^CO 7=10-9 broad component. This ratio is 1.5±0.4 
for the one low-mass YSO for which a decomposition of the 
line profile can be performed in both transitions simultaneously 
(SerSMMl). Thus, it appears that the '^CO 7=10-9 profile be- 
comes increasingly broader compared to the '^CO 7=10-9 pro- 
file with increasing protostellar mass. The average ratio of the 
width of the broad component of the '^CO 7=10-9 line divided 
by the width of the broad component of the '^CO 7=3-2 line 
is approximately 1.0±0.1 for the intermediate-mass sources and 
1.3±0.2 for the low-mass protostars. 

In order to compare the broad velocity component of the 
^^CO data with the narrow C"^0 line profiles across the entire 
studied luminosity range, the FWHM of the '^CO 7=3-2 spec- 
tra is used as a proxy for the FWHM of the '^CO 7=10-9 profiles 
for the high-mass sample. The widths of the fits obtained for the 
^^CO broad velocity components and the narrow C'^O 7=9-8 
and 7=3-2 lines are plotted versus their bolometric luminosi- 
ties (Figs.|4]and|5]l. From the figure of the broad velocity com- 
ponent of the '^CO data we infer that the line-wings become 
broader from low- to high-mass. Low-mass Class protostars 
characterised by powerful outflows, such as L1448, BHR 71 and 
L 1157, are the clear outstanding sources in the plot. The me- 
dian FWHM of this component for each sub-group of protostars 
together with the calculated median of the FWHM values for 
the C"*0 7=3-2 and 7=9-8 Unes are summarised in Table |3] 
Even though there are only six intermediate-mass sources and 
the results could be sample biased, the trend of increasing width 
from low- to high-mass is consistent with the result obtained for 
intermediate-mass objects. 

Regarding the C'^O lines. Fig. |5] and Table [3] show a simi- 
lar behaviour to that observed for the broad component of the 
'^CO but with less dispersion, i.e., the profiles become broader 
from low- to high-mass. This trend is statistically stronger for the 
7=3-2 transition (the Pearson correlation coefficient is higher 
than that calculated for the 7=9-8 line widths) since the number 
of detections is higher for the low-mass sample. The C'^O 7=3- 
2 spectra show slightly narrower profiles than the 7=9-8 line 
for the low- and intermediate-mass sources, with median val- 
ues approximately half the values obtained for the 7=9-8 line 
(see Table [3]). On the other hand, for the high-mass sources the 
median values are practically the same, and similar widths are 
measured for the high-mass sub-sample in both transitions. This 
result will be discussed further in Section 



3.2. Correlations with boiometric iuminosity 

The analysis and characterisation of the line profiles continue 
with the calculation of the integrated intensity of the emission 



Table 3. Median values of the width of the broad velocity com- 
ponent of '^CO 7=10-9 and 7=3-2 spectra, and of the FWHM 
for the C'**0 7=3-2 and 7=9-8 line profiles. 
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Fig. 4. FWHM of the broad velocity component identified in the 
'^CO 7=10-9 line profiles for each type of YSOs versus their 
bolometric luminosities. Low-mass Class (LMO) sources are 
indicated with blue pluses, low-mass Class I (LMI) with black 
triangles, intermediate-mass (IM) YSOs with green asterisks and 
high-mass (HM) objects with red crosses. For the high-mass 
sources, the '^CO 7=3-2 width is used instead. The black dashed 
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line, W= J TiviBdy. This parameter is obtained by integrating the 
intensity of each detected emission line over a velocity range 
which is defined using a 3cr rms cut. 

In order to obtain a more accurate value of W for data with 
lower SIN, such as for the high-7 C'^O lines from the low-mass 
sources, this parameter was approximated to the area of the fitted 
single Gaussian profile. The calculated integrated intensities of 
some sources were compared with measurements from previous 
independent sfiidie s. In the case of NGC 1333 IRAS2A/4A/4B 
dYildizet all 12010). differences in W are not larger than 10%. 
The obtained values from all lines are given in Tables lA.ll to 

Ia3] 

If the emission is optically thin, W is proportional to the col- 
umn density of the specific upper level. In local thermal equilib- 
rium (LTE), the variation of W with 7u characterises the distri- 
bution of the observed sp ecies over the differ ent rotational levels 
(see equation 15.28 from lWilson et al.ll2009l) . In the case of the 
optically thin C'^O 7=9-8 line, the total C'^O column density, 
A^t, is calculated for all sources in order to obtain the H2 col- 
umn density, A^h,- The assu med excitation temp erature, Tex, is 

which shows 



75 K based on the work of lYildiz eta 



io n temp e 
U (1201 oh . 
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Fig. 5. Same as Fig. |4]but for the narrow C'^O 7=9-8 line pro- 
files (top) and C"^0 7=3-2 (bottom). Note that only a few low- 
mass YSOs have been detected in 7=9-8. 



that 90% of the emission in the 7=9-8 transition originates at 
temperatures between 70 and 100 K. The column density Nw^ 
is then obtained by assuming an C'^0/H2 abundance ratio of 
5x10"^. This ratio is obtained by combining the '^O/'^ O iso - 
topologue abundance ratio equal to 540 (Wilson & Roodi l 19941) . 
and the '2CO/H2 ratio as 2.7x10""* (.Lacv et al.,.199 4). The cal- 
culated A^H, values for C'^O 7=9-8 are presented in Table IA.4I 
The integrated intensity is converted to line luminosity, Leo, 
in order to compare these results for sources over a wide rang e 
of distances. The CO and isotopologue line luminosit i es for e ach 
YSO is calculated using equation 2 from IWu et al.l (l2005h as- 
suming a Gaussian beam and point source objects. If the emis- 
sion would cover the entire beam, the line luminosities would 
increase by a factor of 2. The logarithm of this line lumi- 
nosity, log (Leo), is plotted versus the logarithm of the bolo- 
metric luminosity, log(Lboi), for '^CO 7=10-9, ^^CO 7=10- 
9 and C'**0 7=9-8 in Fig. |6l The errors are calculated from 
the rms of the spectrum and considering ~20% distance uncer- 
tainty. A strong correlation is measured (Pearson correlation co- 
efficient r >0.92) between the logarithms of Leo and Lboi for 
all observed CO lines. The top plot of Fig. |6] shows Leo for 
'^CO 7=10-9 emission for all the observed sources versus their 
Lboi. Only one high-mass source was observed as part of WISH 
in this line with HIFI (W3-IRS5) with the value of the inte- 
rated intensity fo r AFGL259 1 obtained from lvan der Wiel et alj 
2013. submittedh . Even though this plot is mainly restricted to 
low- and intermediate-mass sources, a strong correlation is still 
detected over 5 orders of magnitude in both axes. Both low- 
mass Class and Class I YSOs follow the same correlation. 



though the uncertainties of the calculated Leo for these sources 
are higher than for the other types of protostars because the S IN 
is lower All high-mass objects were observed in '^CO 7=10-9, 
so the correlation between log (Leo) and log (Lboi) (Fig- IS] mid- 
dle) is confirmed and extends over almost 6 orders of magnitude 
in both axes. This correlation is also seen for C"*0 7=9-8 but 
with higher dispersion (Fig.|6] bottom), and in the other observed 
transitions of this isotopologue. 

The values of the correlation coefficient and the fit parame- 
ters for all these molecular transitions are presented in Table |4] 
The correlation prevails for all transitions even if the values 
of integrated intensity are not converted to line luminosity. 
Therefore, log(W) still correlates with log (Lboi) over at least 
3 and 6 orders of magnitude on the y and x axis, respectively. 
Similar correlations are obtained when plotting the logarithm 
of Leo versus the logarithm of the source envelope mass, Menv, 
for all the targeted lines (see Fig. |7]for an example using the 
C'**0 7=9-8 line). In these representations, the modelled enve- 
lope mass of the source is directly compared to Leo, a tracer of 
the warm envelope mass. Therefore, the correlation is extended 
and probed by another proxy of the mass of the protostellar sys- 
tem. 

Since the index of the fitted power-law exponents is ~1 
within the uncertainty of the fits (see Table HJi, these correla- 
tions show that log (Leo) is proportional to log (Lboi). In the op- 
tically thin case, this correlation implies that the column den- 
sity of warm CO increases proportionally with the mass of the 
young stellar object. This result can be applied to C'^O because 
the emission lines of this isotopologue are expected to be opti- 
cally thin. Assuming that '^CO 7=10-9 is optically thin as well, 
the column density would increase proportionally with the lu- 
minosity of the source, and practically with the same factor as 
the studied C"*0 transitions. Therefore, even though the condi- 
tions in low-, intermediate- and high-mass star-forming regions 
are different and distinct physical and chemical processes are 
expected to be more significant in each scenario (e.g., ionising 
radiation, clustering, etc.), the column density of CO seems to 
depend on the luminosity of the central protostar alone, showing 
a self-similar behaviour from low- to high-mass. 

To test the optically thin assumption for ^^CO 7=10-9 
and especially for ^■'CO 7=10-9, the line luminosities for the 
'^CO 7=10-9 and t he C^^O 7=10 -9 data were multiplied by a 
i^C/'^C ratio of 65 (Vladilo et al.'T993B and by a i'^0/"*0 ra- 
tio of 540 (Wilson & Rood 1994) , respectively. Therefore, the 
observed and predicted values of Leo for '^CO 7=10-9 and 
'^CO 7=10-9, together with those of C^^^O 7=10-9 can be 
compared across the studied luminosity range (see Fig. |8)- In 
the case of the '^^CO 7=10-9 line, the values of the observed 
and predicted line luminosity are similar (<20% in most of 
them), especially at lower luminosities. In addition, the slope 
of their fits are practically the same within the uncertainty, 
so a similar behaviour is proved. From these results we can 
assume that in general '^CO 7=10-9 is optically thin. For 
'^CO 7=10-9, the ratio of predicted to observed line luminosity 
(65 X Leo[^^CO 7=10-9]/Leo['^CO 7=10-9]) ranges from 0.8 
to 12.5, for IRAS 15398 and W3-IRS5 respectively, and the av- 
erage obtained is 3.3. Therefore, '^CO 7=10-9 is optically thick, 
at least at the line centre which dominates the intensity, and the 
relative value of the optical depth, t, increases slightly with the 
mass of the protostar (r ~1.5 for the low-mass sources, 2.0 for 
the intermediate and ~3.4 for the high-mass object). This is in 



' The '-C/'^C ratio varies with galactocentric radius by up to a factor 
of 2 but this effect is minor and is ignored. 
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(middle) and C'^O 7=9-8 (bottom) spectra for low-mass Class 
(LMO; blue pluses), low-mass Class I (LMI; black triangles), 
intermediate-mass (IM; green asterisks) and high-mass (HM; 
red crosses) YSOs versus their bolometric luminosity. The black 
dashed Une represents the linear function that fits the logarithm 
of the plotted quantities. 



keeping with the expectation that massive YSOs form in the 
densest parts of the giant molecular clouds, GMCs. 

Correcting for optical depth, Lco['^CO 7=10-9] can be used 
to derive Lco['^CO 7=10-9] because both species present a sim- 
ilar behaviour across the luminosity range (similar slopes in their 
fits). This relation can be used in the calculation of Leo for those 
sources for which there are no '^CO 7=10-9 observations, that 
is, the high-mass sample. As highlighted before, using '^CO as 
a proxy for '^CO is restricted to comparisons of integrated in- 
tensities of the emission lines across the studied mass spectrum, 
and cannot be extended to the analysis of the line profile of ^^CO 
and i^CO 7=10-9. 
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Fig. 7. C'^O 7=9-8 line luminosity for low-mass Class 
(LMO; blue pluses), low-mass Class I (LMI; black triangles), 
intermediate-mass (IM; green asterisks) and high-mass (HM; 
red crosses) YSOs versus their envelope masses, M^w- The dash 
black line represents the linear function that fits the logarithm of 
the plotted quantities. 
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Fig. 8. Line luminosity of the '^CO 7=10-9 emission lines, 
red crosses, versus their bolometric luminosity, together with 
the line luminosity of the '^^CO 7=10-9 spectra, blue dia- 
monds, multiply by the assumed abundance ratio of '^C/'^'C for 
the entire WISH sample of YSOs. The line luminosity of the 
C'**0 7=10-9 lines, green triangles, multiply by the assumed 
abundance ratio of '*0/"^0 is plotted together with the pre- 
vious values. The dashed line represents the linear fit of the 
'^CO 7=10-9 spectra, the full line that for the '^^CO 7=10- 
9 transition and the dash-dot line indicates the fit for the 
C^**0 7=10-9 data. 



3.3. Kinetic temperature 

The ratio of the '^CO 7=10-9 and J -3-2 line wings can be used 
to constrain the kinetic temperature Tkin of the entrained outflow 
gas if the two lines originate from the same gas. Yildiz et al. 
(2012) and (2013, submitted) have determined this for the sam- 
ple of low-mass YSOs. Here we consider two sources to inves- 
tigate whether the conditions in the outflowing gas change with 
increasing YSO mass: the intermediate-mass YSO NGC2071 
and the high-mass object W3-IRS5. The critical densities, n„, of 



9 



I. San Jose-Garcfa et al.: HIFI CO observations of YSOs: from low- to high-mass 



Table 4. Slope, b, and intercept, a, of the calculated power-law 
fit for each CO and isotopologue line versus bolometric lumi- 
nosity, together with their errors and the Pearson correlation co- 
efficient, r. The fitted equation is: log (Leo) - a + b- log (Lboi)- 



Line 


a 






h 


r 


'^CO J=\Q-9 


-2.9 ± 


0.2 


0.84 


±0.06 


0.92 


"CO 7=10-9 


-4.4 + 


0.2 


0.97 


± 0.03 


0.98 


C'^O 7=3-2 


-4.1 ± 


0.1 


0.93 


±0.03 


0.98 


C'^O 7=5-4 


-3.5 ± 


0.2 


0.78 


± 0.08 


0.93 


C'^O 7=9-8 


-5.2 ± 


0.2 


1.03 


± 0.05 


0.97 


C'^O 7=10-9 


-5.2 ± 


0.3 


0.96 


± 0.06 


0.96 



the i^CO 7=3-2 and 7=10-9 ti-ansitions at 70 K are -2.0x10^ 
and 4.2x10^ cm"^, respectively. The values were calculated us- 
ing equation 2 from Yang et al. (2010), the CO rate coefficients 
presented in their paper and considering only para-H2 collisions. 
The densities inside the HIFI beam for '^CO 7=10-9 (20") of 
both sources are higher than ria- Therefore, the emission is ther- 
malised and Tkin can be directly constrained by the '^CO 7=10- 
9/7=3-2 line wing ratios. 

The observed ratios of the red and blue wings for these two 
sources as a function of absolute offset from the source ve- 
locity are presented in Fig. [9l These ratios are comp ared with 
the values calculated bv lYildiz et alJ (1201 3. submittedh from the 
'^CO 7=10-9 and 7=3-2 averaged spectra for the low-mass 
Class sample (shaded regions). Since the emission is optically 
thin and we can assume LTE, the kinetic temperatures are cal- 
culated from the equati on ffiat relates the column density and 
the integrated intensity dWilson e t al. 2009). The obtained T^a 
varies from 100 to 210 K. Both the observed line ratios as well 
as the inferred kinetic temperatures are similar to those found 
for the low-mass YSOs, where r^n ranges from 70-200 K for 
Class 0. Although only a couple of higher-mass sources have 
been investigated, the temperatures in the entrained outflow gas 
seem to be similar across the mass range. Note that if part of 
the'^CO 7=10-9 emission originates from a separate warmer 
component, the above values should be regarded as upper limits. 
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Fig. 9. Line-wing ratios of '^CO 7=10-9/7=3-2 as a func- 
tion of absolute offset from the source velocity for the sources 
NGC2071 (top) and W3-IRS5 (bottom) for the red and blue 
wings. The shaded regions are obtained from the a verage spec- 
tra of the low-ma ss sample for these transitions (see lYildiz et al.l 
12013. submittedl) . The spectra have been resampled to 0.6 km s"' 
bin and shifted to km s ' . 



4. Discussion 

The HIFI data show a variety of line profiles with spectra that 
can be decomposed into two different velocity components, such 
as the '^CO 7=10-9 lines, and spectra that show narrow single 
Gaussian profiles (C'^O data). In addition, a strong correlation 
is found between the line and the bolometric luminosity for all 
lines. The next step is to analyse these results to better under- 
stand which physical processes are taking place. 

Section \4~T\ compares the narrow C'^O lines and the '^CO 
broad velocity component in order to better understand the 
physics that these components are tracing and the regions of the 
protostellar environment they are probing. The dynamics of the 
inner envelope-outflow system is studied in Section|42] Finally, 
the interpretation of CO as a dense gas tracer is discussed in 
Sectionfe] 

4. 1 . Broad and narrow velocity components 

The broad velocity component identified in most of the 
'^CO 7=3-2 and 7=10-9 spectra is related to the velocity of 
ffie entrained outflowing material, so th e wings of '^CO can be 
used as tracers of the outflow properties (ICabrit & Bertoutll992t 



i RachiUer & Tafallalll999h . However, there are different effects 
that should be taken into account when this profile component is 
analysed, such as the viewing angle of the protostar and the S/N. 
The former could alter the width of the broad component due to 
projection or even make it disappear if the outflow is located in 
the plane of the sky. Low S/N could also hide the broad compo- 
nent for sources with weak emission. Moreover, the broad veloc- 
ity component should be weaker if the emission lines come from 
sources at later evolutionary stages since th eir outflows become 
weaker and less collimated (see reviews of lBachiller&Tafallal 
1999; Richer et al. 2000; A rce et al. 20'07l) . 

In order to compare ffie line profiles of all observed CO 
lines for each type of YSO and avoid the effects of inclination 
and observational noise playing a role in the global analysis of 
the data, an average spectrum of each line for each sub-type of 
protostar has been calculated and presented in Figs. [TO]and[m 
Regarding the low-mass sample, we observe a striking decrease 
in the width of the broad component from Class to Class I. This 
result shows the decrease of the outflow force for more evolved 
sou rces in ffie low-raass s ample is reflected in the average spec- 
tra (iBontemps et al.lll996h . 
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Fig. 10. i^CO /=10-9 (left) and '^CO 7=10-9 (right) spectra of 
low-mass class (LMO), low-mass Class I (LMl), intermediate- 
mass (IM) and high-mass YSO (HM) averaged independently 
and compared. All spectra were shifted to kms re -binned 
to 0.27 kms"' and the intensity of the emission line scaled 
to a common distance of 1 kpc before averaging. The green 
line indicates the continuum level and the red dashed line the 
kms"' value. W43-MM1 was not included in the average of 
'^CO 7=10-9 high-mass spectra because of the strong absorp- 
tion features caused by HaO^ (see Section l2!2b . 




Fig. 11. Same as Fig.[TO]but for the C'^O spectra from the ob- 
served transitions: 7=5-4 (left), 7=9-8 (centre) and 7=10-9 
(right). More details about these transitions in Appendix lA] 



Higher rotational transitions trace material at higher tem- 
peratures and probe deeper and denser parts of the inner enve- 
lope. For the low- and intermediate-mass sources, the FWHM 
of the C'**0 7=3-2 spectra are generally half that obtained for 
the C"^0 7=9-8 and slightly smaller than those obtained for the 
7=5-4 transition. However, for the high-mass YSOs the values 
of the FWHM are similar for the 7=3-2 and 7=9-8 transitions. 
In order to understand which kind of processes (thermal or non- 
thermal) dominate in the inner regions of the protostellar enve- 
lope traced by our observations, the contribution of these two 
processes to the line width is calculated. The aim is to explain 
whether the broadening of high-7 emission lines is caused by 
thermal or non-thermal motions. 

In the case of the 7=3-2 lines, the upper energy level is 3 1 K, 
so the thermal line width, A^th, is 0.12 kms"' for C'^O at this 
temperature. Comparing this value with the measured FWHM of 
the C'**0 7=3-2 spectra in Table IBTTl we conclude that thermal 
motions contribute less than 5% to the total observed line width, 
Avobs- Therefore, the line width is dominated by non-thermal 
motions Aii„oth- The C'^O 7=9-8 line profiles trace warmer gas 
(up to 300 K) with respect to 7=3-2 increasing the thermal con- 
tribution. However, even at 300 K Auth is 0.68 kms"', which 
means that Aynoth/Auobs is larger than 0.93 even for the low-mass 
sources. Thus, non-thermal motions predominate over thermal 
ones in the studied regions of the protostellar envelopes. These 
motions are assumed to be independent of s cale and do not fol - 
low the traditional size-line width relation (iPineda et al.ll201(ih . 
Therefore, these results are not biased by the distance of the 
source. 

This analysis can be compared to pre-stellar cores, in which 
the line profiles are closer to being dominated by thermal mo- 
tions. For this purpose, the lin e width values cal culated for our 
data are compared to those of iJiiina et al.l (1 1999 1). In that work, 
a database of 264 dense cores mapped in the ammonia lines 
(7, A')=(l,l) and (2,2) is presented. Histograms in Fig. [T2]com- 
pare the values of the line widths observed for pre-stellar cores 
with the observed line width of the C'^O 7=3-2 and 7=9-8 data 
for the WISH sample of protostars for low- (top) and for high- 
mass (bottom) objects. We observe that also for pre-stellar cores, 
the line widths are larger for more massive objects. From these 
histograms and following the previous discussion, we conclude 
that the broadening of the line profile from pre-stellar cores to 
protostars is due to non-thermal motions rather than thermal in- 
crease. Therefore, non-thermal processes (turbulence or infall 
motions) are crucial during the evolution of these objects and 
these motions increase with mass. 



A narrow velocity component has been defined as a line pro- 
file that can be fitted by a Gaussian function with a FWHM 
smaller than 7.5 kms"' (see Section l24l for more details). Since 
C"*0 lines are expected to trace dense quiescent envelope ma- 
terial, high-7 transitions are probing the warm gas in the inner 
envelope. The average C'^O spectra for each type of YSO are 
compared in Fig.[TT] The FWHM of the emission lines increases 
from low- to high-mass protostars (see Fig. |5]and Table |3}. An 
explanation for this result could be that for massive regions, the 
UV radiation from the forming OB star ionises the gas, creating 
an Hii region inside the envelope which increases the pressure 
on its outer envelope. This process may lea d to an increase in the 
turbulent velocity of the envelope material (lMatznerll2002l) . thus 
broadening the narrow component. Therefore, our spectra are 
consistent with the idea that in general, turbulence in the proto- 
stellar envelopes of high-mas s objects is expected to be stronger 
than for low-mass YSOs (e.g. lHeroin et"ani2012h . 



4.2. CO and dynamics: turbulence versus outflow 

'^CO and C'**0 spectra trace diff erent physical stru ctures origi- 
nating close to the protostar (e.g. lYildiz et al.ll2012l) . The broad 
wings of the '^CO 7=10-9 and 7=3-2 data are optically thin 
and trace fast-moving gas, that is, emission from entrained out- 
flow material. On the other hand, the narrow C "*0 spectra probe 
the turbulent and infalling material in the protostellar envelope. 
The relationship between these two different components is still 
poorly understood. 

Following the discussion in § 14.11 we compare the FWHM of 
the narrow component as traced by C '**0 7=9-8 with the FWHM 
of the broad velocity component as traced by '^CO 7=10-9 or 
7=3-2 for the sources detected in both transitions (see Fig. [T3T l. 
The C'^O 7=9-8 data were chosen for this comparison because 
this transition has been observed for the entire sample of YSOs, 
in contrast to C'^'O 7=5^ and 7=10-9 transitions. A correla- 
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Fig. 12. Comparison of the observed width of the ammonia 
line for a s ample of low-mass pre-stellar cores collected by 
iJijina et all d 19991) (green), width of the C"*0 7=3-2 spectra 
(blue) and C'^O /=9-8 line widths (red) for the detected WISH 
sample of low-mass protostars (top). The histogram of the bot- 
tom presents the same values but for high-mass pre-stellar cores 
and the WISH high-mass sample of YSOs. 



tion is found (with a Pearson correlation coefficient r >0.6), in- 
dicating a relationship between the fast outflowing gas and the 
quiescent envelope material. 

Considering the scenario in which the non-thermal compo- 
nent is dominated by turbulence, the relation presented in Fig.[T3] 
indicates that the increase of the velocity of the outflowing gas 
corresponds to an increase in the turbulence in the envelope ma- 
terial and this relation holds across the entire luminosity range. 
One option is that stronger outflows are injecting larger-scale 
movements into the envelope, which increases the turbulence. 
This effect is reflected as a broadening of the C'**0 line width. In 
addition, for the low-mass sources the width of the C'**0 J-9- 
8 lines is larger than for the 7=3-2 spectra, indicating that the 
hotter inner regions of the envelope are more turbulent than 
its cooler outer parts. In the case of the high-mass object, the 
FWHMs of the C'^O 7=9-8 and 7=3-2 lines are comparable, 
which could be partly caused because more luminous YSOs tend 
to be c reated in more m assive and more turbulent molecular 
clouds dWang et al.ll2'009l) . 

Alternatively, infall processes could broaden the C'^0 pro- 
files. Indeed, an increase in FWHM by at least 50% is found for 
C"^0 7=9-8 compared with 7=3-2 in collapsing envelope simu- 
lations due to the higher infall velocities in the inner warm enve- 
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Fig. 13. Width of the C"^0 7=9-8 line profile versus the 
width of the broader component of '^CO 7= 10-9 emission lines 
(crosses) and for the '^CO 7=3-2 line profiles (diamonds). The 
black dashed line represents the linear function that fits the log- 
arithm of the plotted quantities. 



lope dHarsono et al.ll2013. submittedh . If the non-thermal com- 
ponent was dominated by these movements, the sources with 
larger infall rates should have broader C"^0 line widths. These 
objects gene rally have larger outflow rates, i.e., higher o utflow 
activity, (see lTomisakalll998t iBehrend & Maedeill200lh which 
shows up as a broadening of the '^CO wings since the amount 
of material injected into the outflow is larger Therefore, we will 
observe the same result as in the previous scenario, that is, a 
broadening of the '^CO line wings for sources with stronger in- 
falling motions. This relation would hold across the studied lu- 
minosity range. 

Theoretically, large differences in the dynamics of the 
outflow-envelope system between low- and high-mass are ex- 
pected since the same physical processes are not necessarily 
playing the same relevant roles in these different types of YSOs. 
However, Fig. [13] shows that the dynamics of the outflow and 
envelope are equally linked for the studied sample of protostars. 
With the current analysis we cannot disentangle which is the 
dominant motion in the non-thermal component of the line pro- 
file and thus, we cannot conclude if the infalling process makes 
the outflow stronger or if the outflow drives turbulence back into 
the envelope or a combined effect is at play. 

4.3. High- J CO as a dense gas tracer 

The strictly linear relationships between CO high-7 line lumi- 
nosity and bolometric luminosity presented in Figs. |6] to |8] re- 
quire further discussion. The bolometric luminosity of embed- 
ded protostars is thought to be powered by accretion onto the 
growing star and is thus a measure of the mass accretion rate. 
A well-known relation of bolometric luminosity with the out- 
flow momentum flux as measured from '^COJow-7 maps has 
been found across the stellar mass range (iLada&Fichl[l995t 
iBontemps et aP Il996t iRTcher et al.1 l2000l) . so one natural ex- 
planation for our observed '^CO 7=10-9 correlation is that it 
reflects this same relation. However, the outflow wings con- 
tain only a fraction of the '^CO 7=10-9 emission, with the 
broad/narrow intens ity ratios varying from sour ce to source 
(Table lAJ] see also lYildiz et al.ll2013. submittedl for the low- 
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mass sources). Together with the fact that the relations hold 
equally well for '^CO and C'^O high-7 lines, this suggests the 
presence of another underlying relation. Given the strong corre- 
lation with Menv (Fig.|7]l, the most likely explanation is that the 
high-7 CO lines of all isotopologues trace primarily the amount 
of dense gas associated with the YSOs. 

Can this relation be extended to larger scales than those 
probed here? Wu et al. (2005, 2010) found a linear relation be- 
tween HCN integrated intensity and far-infrared luminosity for 
a set of galactic high-mass star forming regions on similar spa- 
tial scales as probed by our data. They have extended this re- 
lation to include extragalactic sources to show that this linear 
regime ex tends to the scale s of en tire galaxies, as first demon- 
strated by iGao & SolomonI (|2004 . In contrast, the CO 7=1-0 
line shows a superlinear relation with far-infrared luminosity 
(sometimes converted into star formation rate) and with the total 
(H I + H2) gas su rface density with a power-law index of 1.4 
(see lKennicutt & Evans.2012 for review). The data presented in 
this paper suggest that it is the mass of dense molecular gas as 
traced by HCN that controls the relation on large scales rather 
than the mass traced by CO 7=1-0 (see also Lada et al. 2012). 
With the increased number of detections of '^CO 7=10- 
local and high re dshift g alaxies with Herschel 



y m local and high re dshiit g alaxies with He rschei 
(e.g.. lyaiTderWerfet aLl |2010t ISpinoglio et ail 120121: 
iKamenetzkY et al. 2012 *: 'Meiierink et al."2013') and millimeter 



interferometers (e.g., fWeiB et al. 2007; Scott et al. 2Q\%, the 
question arises whether this line can serve as an equally good 
tracer of dense gas. As an initial test, we present in Fig. [14] our 
'^CO 7=10-9 - LpiR (assuming Lboi = T-fir) relation with these 
recent extragalactic detections included. As can be seen, the 
relation does indeed extend to larger scales. Given the small 
number statistics, the correlation has only limited meaning and 
this relation needs to be confirmed by additional data. 

An alternative v iew has been presente d b 

iKrumholz & Thompson! (l2007h and iNaravanan elal] (I2OO 
who argue that the linearity for dense gas tracers is a coinci- 
dence resulting from the fact that these higher excitation lines 
are subthermally excited and probe only a small fraction of the 
total gas. They suggest that the relation should change with 
higher critical density tracers and even become sublinear with a 
power-law index less than 0.5 for transitions higher than 7=7-6 
for CO. On the small scales probed by our data, we do not see 
this effect and the linear relations clearly continue to hold up to 
the 7=10-9 transition. 



5. Conclusions 

The analysis of the ^^CQ 7=3-2, 7=10-9, '^CO 7=10-9, 
C"^0 7=3-2, 7=5-4, 7=9-8 and 7=10-9 line profiles allows 
us to study several fundamental parameters of the emission line, 
such as the line width and the line luminosity, and to constrain 
the dynamics of different physical structures of WISH protostars 
across a wide range of luminosities. Complementing the HIFI 
data, lower-7 observations from the JCMT are included in order 
to achieve a uniform picture of the interaction of YSOs with their 
immediate surroundings. Our results are summarised as follows: 

• A gallery of line profiles identified in the HIFI CO spectra 
is presented. '^CO and '^CO 7=10-9 line profiles can be 
decomposed into two different velocity components, where 
the broader component is thought to trace the entrained out- 
flowing gas. This broad component weakens from the low- 
mass Class to Class I stage. Meanwhile, the narrow C'^O 
lines probe the bulk of the quiescent envelope material. The 
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Fig. 14. Line luminosity of the '^CO 7=10-9 emission versus 
their bolometric luminosity for the entire WISH sample of YSOs 
(same notation as Fig. |6]l and for four galaxies in our local uni- 
verse and two at high redshifts (black diamonds). The solid line 
represents the linear fit of the galactic '^CO 7=10-9 spectra and 
the dashed line the linear fit including the values of the six galax- 
ies. 

widths are dominated by non-thermal motions including tur- 
bulence and infall. The next step will be to constrain the con- 
tribution of these non-thermal mechanisms on the C'**0 line 
profiles by using radiative transfer codes such as RATRAN 
(Hogerheiide & van der Tak 2000). 

• The narrow C'**0 7=9-8 line widths increase from low- to 
high-mass YSOs. Moreover, for low- and intermediate-mass 
protostars, they are about twice the width of the C'^O 7=3- 
2 lines, suggesting increased turbulence/faster infall in the 
warmer inner envelope compared to the cooler outer enve- 
lope. For high-mass objects the widths of the 7=9-8 and 
7=3-2 lines are comparable, suggesting that the molecular 
clouds in which these luminous YSOs form are more turbu- 
lent. Extending the line width analysis to pre-stellar cores, a 
broadening of the line profile is observed from these objects 
to protostars caused by non-thermal processes. 

• A correlation is found between the width of the '^CO 7=3-2 
and 7=10-9 broad velocity component and the width of the 
C"^0 7=9-8 profile. This suggests a link between entrained 
outflowing gas and envelope motions (turbulence and/or in- 
falling) which holds from low- to high-mass YSOs. This 
means that the interaction and effect of outflowing gas and 
envelope material is the same across the studied luminos- 
ity range, indicative of the existence of an underlying com- 
mon physical mechanism which is independent of the source 
mass. 

• A strong linear correlation is found between the logarithm 
of the line and bolometric luminosities across six orders of 
magnitude on both axes, for all lines and isotopologues. This 
correlation is also found between the logarithm of the line 
luminosity and envelope mass. This indicates that high-7 CO 
transitions (up to 7=10-9) can be used as a dense gas tracer, 
a relation that can be extended to larger scales (local and high 
redshift galaxies). 
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Appendix A: Characterisation of the IHIFI data 

The main characteristics of the HIFI data, together with the 
spectra are presented in this appendix. The description of the 
observed lines focuses first on '^CO 7=10-9 (Section lA.ll ). 
where the main characteristics for the low-, intermediate- and 
high-mass sources are listed in this order Next we discuss the 
^^CO 7=10-9 spectra (Section [A. 2l i. following the same struc- 
ture, and finally all observed C'**0 lines (Section [A. 3l i are pre- 
sented. 

A.1. ^^CO J=10-9 line profiles 

The '^CO 7=10-9 line was observed for the entire sample of 
low- and intermediate-mass YSOs and for one high-mass ob- 
ject (W3-IRS5). All the observed sources were detected and the 
emission profiles are the strongest and broadest among the tar- 
geted HIFI CO lines (see Fig. lA. II and Table lA.ll for further in- 
formation). Within the low-mass Class sample, the main beam 
peak temperature, T^*, ranges from 0.8 to 8.1 K. As indicated 
in Section 1231 the emission in some sources (around 73 % of 
the detected lines of this sub-group) can be decomposed into 
two velocity components. The FWHM of the narrower compo- 
nent varies from 2.3 kms"' to 9.3 kms"', while the width of 
the broad component shows a larger variation, from 8.3 kms"' 
to 41.0 kms '. For the low-mass Class I YSOs, similar inten- 
sity ranges are found, with varying from ~0.6 to 10.2 K. 
The Class I protostars present narrower emission lines than the 
Class sources and only the 27 % of the emission line profiles 
can be decomposed into two velocity components. The narrow 
component ranges from 1.8 to 5.1 kms ' and the broad compo- 
nent varies from 10.0 to 16.1 kms 

For the intermediate-mass protostars, the intensity increases 
and varies from T^^=2A to ~28.0 K. The profiles are broader 
and the distinction between the two different components is 
clearer than for the low-mass sources, so all profiles can be fitted 
with two Gaussian functions. The FWHM of the two identified 
velocity components varies from 2.7 to 7.6 km s"' for the narrow 
component and from 15.6 to 24.3 kms"' for the broad compo- 
nent. 

Only W3-IRS5 was observed in '^CO 7=10-9 from the 
high-mass sample. The spectrum is presented in Appendix |B] 
Fig. IB. 31 together with other lines of this source. The profile is 
more intense than any of the low- and intermediate-mass sources 
(^mb''=48.5 K) and has the largest FWHM; 8.4 kms"' for the 
narrower component and 28.6 km s ' for the broad component. 

Self-absorption features have been detected in 5 out of 33 
obsei-ved '^CO 7=10-9 emission lines (the sources are indicated 
in Table lA.ll . However, these features are weak and of the order 
of the rms of the spectrum, so no Gaussian profile has been fit- 
ted. No specific symmetry can be determined, that is, there is no 
systematic shift in the emission of the broad component relative 
to the source velocity (see Fig. lA.ll for comparison). Overall, the 
data do not show any infall signature. 

A.2. CO J =10-9 line profiles 

The observed '^CO 7=10-9 emission lines for the low- and 
intermediate-mass sources are less intense, narrower and have 
a lower S/N than the '^CO 7=10-9 spectra (Fig.|A21l. TablelO 
contains the parameters obtained from the one or two Gaussian 
fit to the detected line profiles. In the case of the low-mass 
Class spectra, three sources are not detected down to 17 mK 



rms in 0.27 kms bins and the profile of only two sources 
(SerSMMl and NGC1333 IRAS4A) can be decomposed into 
two different velocity components. The T"^^ ranges from 0.05 
to 0.8 K and the FWHM of the narrow profiles varies from 0.7 
to 6.8 km s"', with the largest values corresponding to the broad 
velocity component being 13.2 kms"' for NGC1333 IRAS4A. 
In the case of the Class I sample, four sources are not detected 
and none of the emission line profiles can be decomposed into 
two velocity components. The averaged intensity is lower than 
for the Class objects, ranging from 7'^'^g''=0.05 to 0.52 K. The 
value of the line width also drops and the interval varies from 
1.5 to 7.3 kms"'. 

For the intermediate-mass YSOs, a better characterisation 
of the line profile is possible since the lines are stronger and 
have higher S/N than the low-mass objects with ranging 
from 0.1 to 2.7 K. Compared to the '-CO 7=10-9 profiles, 
the '^CO 7=10-9 lines are more symmetric and only the emis- 
sion profile of one source can be decomposed into two different 
Gaussian components. Regarding the FWHM of the lines fitted 
by the narrow Gaussian, the interval goes from 4.3 to 6. 1 km s"' . 

Around 63% of the detected '^CO 7=10-9 emission lines 
(12 out of 19) for the high-mass YSOs can be decomposed 
into two distinct velocity components, whereas the decompo- 
sition of the profiles is only possible for 10% of the detected 
low-mass objects (2 out of 20) and for ~17% of the detected 
intermediate-mass YSOs (1 out of 6). The reason for the lower 
percentage recorded for the low- and intermediate-mass sources 
could be caused by a lower S/N compared with the bright high- 
mass sources. The weakest line from the high-mass sample has 
a TT!'^ of 0.7 K and the most intense a rff^'' of 20.8 K. The 

Md Mrs 

FWHM of the narrower component varies from 3.3 kms"' to 
7.2 kms"'. The width of the broad component presents a larger 
variation since the minimum value is 8.7 kms"' and the maxi- 
mum 21.9 kms"'. This component appears either red- or blue- 
shifted. There is no significant trend with evolution stage as 
probed by the presence of IR-brightness or ionising radiation 
(Fig.|A3. 

A.3. C^^O line profiles 

Three transitions of C'**0 were obtained within WISH together 
with water observations: C'**0 7=5-4, 7=9-8 and 7=10-9. 
Only Class and intermediate-mass YSOs were observed in 
C'**0 7=5-4, tracing regions with an upper energy level of 
~79 K (see Fig. IA.3I ). This line is obtained in parallel with a 
deep integration on the 548 GHz H^^O l\o - loi transition for 
19 sources. Thus, the spectra have very high S/N with an rms 
of 9 mK for low-mass Class sources and less than 20 mK for 
intermediate-mass YSOs in 0.27 kms"' bins. The main char- 
acteristic of this transition is the narrow profile seen in all the 
emission lines for the narrower component, with a FWHM of 
less than 2.0 km s"' for the low-mass sources, and 3.7 km s"' for 
the intermediate-mass objects. In addition, other features are de- 
tected thanks to the high S/N, e.g. a weak broad velo city compo- 
nent for the low-mass objects NGC1333 IRAS 4A ( Yildiz et al.l 
[2010), L483, SerSMMl and SerSMM4. This componentis also 
identified in the C'**0 7=5^ line for the intermediate-mass 
sources NGC2071 (see Fig.O and Vela IRS 19. The values of the 
single or two Gaussian fit of these lines are presented in Table 
IA3] 

For the 7=9-8 transition, ~55% of the observed lines are 
detected, probably due to the lower S /N caused by shorter ex- 
posure times than for the 7=5-4 line. The lines are detected 
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Fig. A.l. '^CO 7=10-9 spectra for low- and intermediate-mass YSOs. The green line represents the baseline level and the pink 
Gaussian profile the broad velocity component for those sources for which a 2 Gaussian decomposition has been performed. All the 
spectra have been shifted to zero velocity. The numbers indicate where the spectra have been scaled for greater visibility. 



in 5 out of 26 low-mass sources; 4 out of 6 intermediate-mass 
YSOs; and in all 19 high-mass protostars (see Fig. IA.4) . The 
C'^O 7=9-8 emission lines appear weak with median val- 
ues of 0.10, 0.14 and 0.83 K for the low-, intermediate- and 
high-mass objects, respectively. Most of the emission line pro- 
files of C'^O 7=9-8 can be fitted by a single Gaussian with 
a FWHM from 2.0 kms"' to 3.9 kms"' for the low-mass ob- 
jects; from 2.8 kms"' to 5.4 kms"' for the intermediate-mass 
sources; and from 3.1 to 6.4 kms ' for the high-mass YSOs 
(values summarised in Table lA.4l i. Only a two Gaussian decom- 
position has been performed for three ultra-compact HII regions 
(G10.47-H0.03,W51N-el and G5.89-0.39). For these objects, the 
broad velocity components are larger than 1 6 km s " ' , as is shown 
inFig.lAH 

Finally, the C'**0 7=10-9 transition was observed in 
30 minute exposures for all low-mass Class protostars, 
one intermediate-mass source and all high-mass YSOs. 



Additional deeper integrations of 300 minutes were ob- 
tained for NGC1333IRAS2A (as part of WISH) and 
for NGC1333IRAS4A, NGC1333 IRAS 4B, Elias 29 and 
GSS30IRS1 (as part of open-time programme OT2_rvisser_2) 
in parallel with deep Hj^O searches. The line was detected in 
five low-mass sources and in all 19 high-mass objects (Fig. IA.5l l. 
This line appears close to the 1097 GHz H2O 3i2-3o3 transition. 
For most of the high-mass objects, the line profile of this water 
transition shows broad wings which extend few kms so the 
C'^O 7=10-9 emission line is found on top of the broad water 
red wing. In order to properly analyse the emission of this CO 
isotopologue, the line wings of the 1097 GHz water transition 
were fitted with a Gaussian profile, subtracted and the residu- 
als plotted. With this method, the C'^O 7=10-9 emission line 
for the high-mass sample has been isolated. The temperature of 
the gas that 7=10-9 traces is likely similar to that traced by the 
7=9-8 transition, so the lines are also weak with median values 
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Table A.l. Observed and fitted properties of the '^CO 7=10-9 line profiles. The subscript B refers to the broad velocity component 
and the to the narrow component. 



Source 


rms" 


/ TuBdv" 


^ MBB 


^ MBN 


. ,peak 


, rpeak 


FWHMb 
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(mK) 


(Kkms-') 
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(km s-') 


(km s-') 


(kms-') 


(km s-') 
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L 1448-MM'' 


91 


21.5 


0.4 


0.9 


10.3 


6.0 


41.0 


4.8 


NGC 1333 1RAS2A 


104 


20.3 


0.4 


1.5 


12.8 


8.2 


8.3 


3.9 


NGC 1333 1RAS4A" 


105 


45.5 


1.8 


- 


10.3 


- 


24.8 


- 


NGC 1333 1RAS4B 


104 


32.4 


1.4 


1.5 


8.1 


7.1 


16.5 


3.3 


L 1527'' 


93 


4.8 


- 


1.6 


- 


4.9 


- 


2.4 


Cedll01RS4'^ 


127 


4.9 


- 


1.6 


- 


4.2 


- 


2.5 


IRAS 15398 


132 


16.5 


0.5 


2.5 


-1.0 


4.1 


15.0 


4.2 


BHR71'''' 


111 


9.9 


0.2 


0.6 


-8.5 


-5.4 


30.8 


9.3 


L483-MM 


108 


10.9 


0.3 


1.5 


2.8 


5.3 


19.1 


2.9 


Ser SMM 1" 


98 


81.3 


3.4 


4.1 


7.1 


8.6 


15.6 


5.9 


Ser SMM 3 
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33.7 


1.0 


2.0 


8.8 


7.0 


20.8 


4.1 


Ser SMM 4 


97 


39.0 


1.6 


3.3 


2.1 


6.8 


10.7 


4.3 


L723-MM'- 


110 


6.8 


- 


1.1 


- 


10.9 


- 


4.9 


B335 


120 


11.1 


0.5 


1.0 


8.9 


8.3 


16.1 


2.3 


LI 157 


103 


9.5 


0.3 


0.5 


0.4 


2.9 


26.8 


2.7 


Low-mass: Class I 
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0.4 
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Vela IRS 19 


110 
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14.9 
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22.0 


3.3 


Vela IRS 17 
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94.0 
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5.6 
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7.0 


13.9 


8.2 


10.9 


24.3 


7.6 


AFGL490 


129 


29.6 


1.2 


3.0 


-11.5 


-13.4 


15.6 


2.7 


High-mass 


















W3-IRS5'' 


102 


674.8 


10.8 


32.0 


-40.7 


-37.5 


28.6 


8.4 



Notes. In 0.27 km s ' bins. Integrated over the entire line, not including "bullet" emission. Single Gaussian fit. '''' EHV emission features 
removed from the spectra by using two additional Gaussian fit profiles. Self-absorption features detected. 



of of 0.04 K for the low-mass protostars and 0.52 K for 
the high-mass objects. The FWHM of the one single Gaussian 
profile which fits these lines are slightly larger, ranging from 3.4 
to 7.5 km s~' for the high-mass sources (Table lA3] for more de- 
tails). 



Appendix B: JCMT data 

The central spectrum of the '^CO and C'^O J=3-2 spectral maps 
observed with the HARP instrument of the JCMT are presented 
in this section (see Figs. lB.4l and lB.5b . These central spectra were 
convolved to a 20" beam in order to compare them with the 
HIFI data. The '^CO and C'^O 7=3-2 spectra for the low-mass 
sources BHR71, CedllOIRS4, IRAS 12496 and HH46 were 
observed with APEX because of their low declination. Similarly 
to the JC MT data, the central spectrum w as convolved to a 20" 
beam. See lYildiz et al.l ( 12013. submittedh for more information 
about the low-mass protostar observations. 



The spectral maps are not presented in this paper because 
only the values of FWHM from the broad component for the 
'^CO 7=3-2 central spectrum were used in the analysis and dis- 
cussion, together with the width and integrated intensity of the 
C'**0 7=3-2 data. These values are presented in Table IB. Il Some 
spectra from these species are also plotted together with the HIFI 
data in the figures shown in this appendix in order to make a di- 
rect comparison of the lines for different types of protostars. 

The '^CO 7=3-2 data show more complex line profiles 
than the 7=10-9, as indicated in Section [TT] with intense self- 
absorption features and broad velocity components (see spectra 
from Fig. IB.4b . More than 82% of the observed and detected 
lines (39 out of 47) present a broad velocity component, and the 
FWHM ranges from 7.4 to 53.3 kms"', values corresponding to 
a low- and high-mass YSOs respectively. 

On the other hand, the narrow C '^O 7=3-2 spectra show sin- 
gle Gaussian emission profiles, similar to those from higher-7 
transitions observed with HIFI. The FWHM of these data varies 
from 0.6 to 7.3 kms The spectra is presented in Fig. IB. 51 
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Table A.2. Observed and fitted properties of the '^'CO /=10-9 line profiles for the detected sources. The subscript B refers to the 
broad velocity component and the to the narrow component. 



Source 


rms" 




rT-.peak 
^ MBB 


ry-Tpeak 
MBN 


. ,peak 


. rpeak 


FWHMb 


FWHMfi 




(mK) 


(Kkms-') 


(K) 


(K) 


(kms-') 


(kms-') 


(km s-') 


(kms-') 


Low-mass: Class 


















L 1448-MM 


20 


0.2 


- 


0.1 


- 


4.7 


- 


1.6 


NGC 1333 IRAS 2 A 


15 


0.8 


- 


0.2 


- 


7.5 


- 


2.1 


NGC 1333 IRAS 4A 


23 


1.1 


0.07 


0.09 


5.7 


6.7 


13.2 


0.7 


NGC 1333 IRAS4B 


17 


0.8 


- 


0.1 


- 


6.9 


- 


6.8 


IRAS 15398 


24 


0.2 


- 


0.16 


- 


5.0 


- 


1.1 


BHR71 


17 


0.4 


- 


0.2 


- 


-4.6 


- 


1.8 


L483-MM 


16 


0.3 


- 


0.1 


- 


5.0 


- 


2.4 


Ser SMM 1 


22 


4.0 


0.2 


0.5 


7.6 


8.3 


10.6 


2.4 


Ser SMM 3 


20 


0.3 


- 


0.1 


- 


7.2 


- 


4.0 


Ser SMM 4 


25 


0.6 


- 


0.1 


- 


5.7 


- 


6.0 


L 723-MM 


16 


0.3 


- 


0.1 


- 


11.5 


- 


2.7 


B335 


20 


0.3 


- 


0.1 


- 


8.1 


- 


2.3 


LI 157 


20 


0.1 


- 


0.05 


- 


3.0 


- 


2.2 


Low-mass: Class I 


















L1489 


32 


0.4 


- 


0.05 


- 


7.8 


- 


7.3 


L1551 IRS 5 


29 


1.3 


- 


0.47 


- 


6.5 


- 


2.4 


TMRl 


29 


0.3 


- 


0.08 


- 


6.0 


- 


3.6 


HH46 


30 


0.1 


- 


0.11 


- 


5.4 


- 


1.5 


IRAS 12496 


28 


0.8 


- 


0.15 


- 


3.1 


- 


4.2 


GSS30IRS1 


33 


2.0 


- 


0.46 


- 


2.8 


- 


3.4 


Elias 29 


28 


2.1 


- 


0.38 


- 


4.7 


- 


5.1 


Intermediate-mass 


















L1641S3MMS1 


38 


1.8 


- 


0.2 


- 


4.7 


- 


6.1 


Vela IRS 19 


41 


1.1 


- 


0.2 


- 


12.1 


- 


5.9 


Vela IRS 17 


43 


4.9 


- 


0.8 


- 


0.4 


- 


4.4 


NGC 7 129 FIRS 2 


20 


0.9 


- 


0.1 


- 


-9.6 


- 


4.3 


NGC 2071 


17 


16.6 


0.4 


2.1 


9.1 


9.8 


14.7 


4.7 


AFGL490 


47 


2.6 


- 


0.3 


- 


-13.5 


- 


4.8 


High-mass 


















IRAS05358-)-3543 


18 


6.2 


0.2 


0.7 


-15.1 


-15.9 


15.4 


3.7 


IRAS 16272-4837 


18 


3.9 


- 


0.7 


- 


-47.0 


- 


4.8 


NGC6334-I(N) 


23 


13.8 


0.5 


1.3 


-5.7 


-3.8 


13.3 


4.4 


W43-MM1 


38 


4.5 


- 


0.8 


- 


98.4 


- 


5.8 


DR21(OH) 


18 


52.2 


0.9 


5.7 


-1.9 


-3.3 


16.0 


5.8 


W3-IRS5 


21 


129.6 


2.1 


17.0 


-38.7 


-38.3 


14.7 


4.7 


IRAS18089-I732 


27 


8.3 


- 


1.4 


- 


33.0 


- 


4.6 


W33A 


19 


11.4 


0.5 


1.1 


36.8 


37.7 


11.7 


3.5 


IRAS18151-I208 


19 


3.4 


- 


0.8 


- 


33.4 


- 


3.6 


AFGL2591 


19 


29.5 


0.8 


5.4 


-5.8 


-5.4 


9.8 


3.3 


G327-0.6 


25 


12.8 




1.8 




-44.8 




6.6 


NGC6334-I 


23 


57.7 


0.7 


7.2 


-6.7 


-6.8 


17.8 


5.4 


G29.96-0.02 


46 


31.3 


0.6 


4.6 


96.9 


98.8 


13.6 


4.6 


G31.4I+0.3I 


52 


17.4 




2.5 




96.8 




6.1 


G5.89-0.39 


34 


133.7 


2.6 


6.2 


10.5 


9.4 


21.9 


5.9 


GIO.47+0.03 


200 


26.2 


1.8 




66.8 




10.0 




G34.26-i-0.15 


59 


58.8 


1.3 


6.7 


58.0 


57.7 


12.0 


5.0 


W5IN-el 


50 


70.2 


1.3 


5.0 


58.6 


57.1 


18.6 


7.2 


NGC7538-IRS1 


24 


41.0 


2.0 


5.3 


-59.8 


-57.6 


8.7 


3.5 



Notes. In 0.27 km s ' bins. Integrated over the entire line. 



and the constrained values of FWHM and integrated intensity 
in Table lO 
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Table A.3. Observed and fitted properties of the C^^O J=5-4 line profiles for the observed sources. The subscript B refers to the 
broad velocity component and the to the narrow component. 



Source 


rms" 




^ MBB 


MBN 


. .peak 


. .peak 
N 


FWHMb 






(mK) 


(Kkms ') 


(K) 


(K) 


(kms^') 


(kms^') 


(kms^') 


(kms^') 


Low-mass: Class 


















L 1448-MM 


4 


0.47 


- 


0.37 


- 


5.1 


- 


1.3 


NGC1333 IRAS2A 


2 


0.87 


- 


0.52 


- 


7.5 


- 


1.4 


NGC1333 IRAS4A 


4 


0.68 


0.03 


0.32 


7.9 


6.9 


7.8 


1.3 


NGC1333 IRAS4B 


3 


0.32 


- 


0.14 


- 


7.0 


- 


1.9 


L1527 


3 


0.55 


- 


0.23 


- 


5.7 


- 


2.0 


BHR71 


3 


0.63 


- 


0.37 


- 


-AJ 


- 


1.4 


L483-MM 


3 


0.31 


0.03 


0.16 


5.2 


5.4 


4.3 


1.1 


Ser SMM 1 


4 


1.93 


0.13 


0.84 


8.7 


8.4 


4.2 


1.4 


Ser SMM 3 


5 


1.40 


- 


0.55 


- 


7.6 


- 


2.0 


Ser SMM 4 


3 


1.55 


0.13 


0.44 


7.6 


7.7 


4.7 


1.8 


L723-MM 


5 


0.19 


- 


0.09 


- 


10.0 


- 


1.9 


B335 


4 


0.26 




0.22 




8.2 




1.2 


L1157 


4 


0.15 




0.09 




2.6 




1.3 


Intermediate-mass 


















L1641 S3MMS1 


20 


0.73 




0.40 




5.2 




1.7 


Vela IRS 19 


10 


1.67 


0.05 


0.46 


12.5 


11.6 


9.2 


2.4 


Vela IRS 17 


25 


5.32 




1.25 




4.1 




3.7 


NGC 7 129 FIRS 2 


4 


0.55 




0.18 




-9.9 




2.3 


NGC2071 


5 


7.29 


0.55 


1.33 


8.9 


9.6 


7.4 


2.0 


AFGL490 


14 


3.55 




1.04 




-13.3 




2.7 



Notes. In 0.27 km s ' bins. Integrated over the entire line. 
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Low— mass: Class 

-'sMMl I 



- IRAS 


2A 









- L1551 


IRS5 - 
xO.5 





- SMM3 


x2 - 







- L115 

JjIL 


? 

AhjIiIi 











- IRAS4I 


3 





- SMM4 










BHR71 



CedllO IRS4 



Low— mass: Class I 



- Oph IRS63 



^ril ■nnril|li|-iii[ 




- GSS30 


IRSl - 
xO.4- 





- TMCl 




fcii ill, 









- HH46 

Mini 


i.Al 









1i,i4A,Ja« aM.l 



IRAS12496 



I 30 -30 

V (km/s) 



30 -30 30 

V (km/s) 



30 -30 

V (km/s) 



Intermediate— mass 




: NGC7129 


FIRS2 







: NGC 2071 







: AFGL490 




: ; 


x4 1 





-20 -10 10 20 -20 -10,, ,10, 20 -2" ->» » m 20 

V 



0, ,10 

(km/ s) 

High— mass 




NGC6334-I(N) 






M 







"W33A 















IRAS18151-1Z08 



"ArGL2591 







;G3Z7-o.e 


\ x3 - 





;NGC6334 


1 -I ' 





^GS.ag- 


Ao.39 ' 









GlO.47 + 0.03 





20 -20 20 -30 30 -20 20 -20 20 

V (km/s) 



Fig.A.2. Same as Fig. lA.ll but for the '^CO 7=10-9 spectra from the observed low-, intermediate- and high-mass YSOs. In the 
figure with the high-mass sample, the sources are presented according to their evolutionary stage. In the first row we find the mid- 
IR-quiet high-mass protostellar objects (HMPOs), in the second row the mid-IR-bright HMPOs, in the third line hot molecular cores 
and in the last row the ultra-compact HII regions. 
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Table A.4. Observed and fitted properties of the C'^O 7=9-8 narrow line profiles for the detected sources. 



Source 


/ TmbcIv* 


^ MB 


FWHM 


rms" 






(K kms ) 


(K) 


(kms"') 


(mK) 


(cm~^) 


Low-mass: Class 












NGC 1333 IRAS2A 


0.21 


0.07 


2.0 


20 


1.2x10^' 


Ser SMM 1 


0.70 


0.13 


3.0 


23 


4.1x10^' 


Low-mass: Class I 












T 1 C C 1 Tn o c 

L 1551 IRS 5 


0.22 


0.10 


3.1 


24 


1.3x10 


GSS30IRS1 


0.23 


0.06 


2.2 


25 


1.4x10^' 


Elias 29 


0.40 


0.08 


3.9 


30 


2.3x10^' 


Intermediate-mass 












Vela IRS 19 


0.37 


0.11 


3.2 


24 


2.2x10^' 


Vela IRS 17 


0.69 


0.16 


2.8 


25 


4.0x10^' 


NGC 2071 


2.97 


0.41 


4.6 


25 


1.7x10^^ 


AFGL490 


0.45 


0.08 


5.4 


20 


2.6x10^' 


High-mass 












IRAS05358+3543 


0.93 


0.22 


4.0 


56 


5.4x10^' 


IRAS 16272-4837 


1.18 


0.14 


5.9 


55 


6.9x10^' 


NGC6334-I(N) 


1.82 


0.37 


3.9 


55 


1.1x10^- 


W43-MM1 


1.81 


0.29 


5.9 


49 


1.0x10^^ 


DR21(0H) 


9.11 


1.51 


5.5 


57 


5.3x10^^ 


W3-IRS5 


25.43 


4.80 


4.2 


66 


1.5x10^' 


IRAS18089-1732 


3.26 


0.52 


4.1 


53 


1.9x10^^ 


W33A 


3.17 


0.43 


4.8 


46 


1.9x10^^ 


IRAS18151-1208 


0.60 


0.24 


4.2 


54 


3.5x10^' 


AFGL2591 


4.56 


1.21 


3.1 


39 


2.7x10^^ 


G327-0.6 


6.16 


0.83 


6.1 


57 


3.6x10^^ 


NGC6334-I 


10.97 


2.00 


5.0 


54 


6.4x10^2 


G29.96-0.02 


7.50 


1.41 


4.3 


41 


4.4x10^^ 


G31.41+0.31 


5.26 


0.73 


6.0 


36 


3.1x10^^ 


G5. 89-0.39'' 


31.35 


2.13 


5.9 


51 


1.8x10^3 


GlO.47+0.03'' 


5.71 


0.53 


6.4 


79 


1.1x10^^ 


G34.26+0.15 


16.26 


2.53 


5.5 


49 


9.5x10^2 


WSlN-cl'' 


15.22 


1.83 


5.7 


48 


8.9x10^- 


NGC7538-1RS1 


6.32 


1.37 


4.0 


48 


3.7x10" 



Notes. In 0.27 km s ' bins. Integrated over the entire line. Column density of H2 calculated for an excitation temperature of 75 K and a 
C'^0/H2 abundance ratio of 5x10"^. Two Gaussian fit. Only the values from the narrow velocity component are presented. 
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Low— mass: Class 




- SMMl 
■ j 


1 





- L723-MM 




- IRAS4 


A 

1 x2 





- SMM3 







- L4B3 










- SMM4 

: J 







-10 10 



1 


- L1448 


MM 


0,5 




x2 - 










-10 10 



-10 10 -10 10 



V (km/s) V (km/s) 

Intermediate — mass 




Table A,5. Observed and fitted properties of the /=10-9 
line profiles for the detected sources. 



Source 


/ TuBdv'' 




FWHM 


rms" 




(Kkms-') 


(K) 


(km s-') 


(mK) 


Low-mass: Class 










NGC 1333 IRAS 2A 


0.34 


0.04 


10.7 


8 


NGC 1333 IRAS 4B 


0.05 


0.03 


1.4 


10 


Ser SMM 1 


0.61 


0.07 


4.8 


17 


Low-mass: Class I 










GSS30IRS1 


0.12 


0.04 


3.1 


8 


Elias29 


0.22 


0.04 


4.7 


7 


High-mass 










IRAS05358+3543 


0.42 


0.10 


3.8 


26 


IRAS 16272-4837 


0.49 


0.11 


4.8 


26 


NGC6334-I(N) 


2.31 


0.28 


7.5 


20 


W43-MM1 


0.92 


0.12 


5.3 


23 


DR21(OH) 


7.94 


1.08 


6.0 


33 


W3-IRS5 


17.04 


3.39 


4.1 


22 


IRAS18089-1732 


1.54 


0.30 


4.5 


28 


W33A 


1.72 


0.30 


4.3 


22 


IRAS18151-1208 


0.48 


0.09 


4.2 


22 


AFGL2591 


3.80 


0.88 


3.4 


23 


G327-0.6 


3.03 


0.42 


5.8 


34 


NGC6334-I 


10.09 


1.55 


5.3 


28 


G29.96-0.02 


5.79 


1.02 


4.5 


30 


G31.41+0.31 


3.98 


0.52 


6.1 


38 


G5.89-0.39 


11.59 


2.08 


7.0 


27 


GlO.47+0.03 


3.44 


0.47 


7.1 


26 


G34.26+0.15 


11.80 


1.73 


6.0 


29 


W51N-el 


10.49 


1.48 


6.8 


52 


NGC7538-IRS1 


3.95 


0.88 


3.9 


37 



Notes. In 0.27 km s ' bins. **' Integrated over the entire line. 



V (km/s) 



Fig. A.3. Same as Fig. lA.ll but for the C'**0 7=5-4 spectra of 
the low-mass (Class 0) and intermediate-mass protostars. For the 
low-mass sample the HRS spectra are presented while for the 
intermediate-mass objects the WBS data is used. 
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Fig.A.4. Same as Fig. lA.ll but for the C'**0 7=9-8 spectra from the low-, intermediate- and high-mass YSOs. The high-mass 
objects are presented according to their evolutionary stage, as is explained in Fig. lA.2l 
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Fig. A.5. Same as Fig. lA.ll but for the C'**0 7=10-9 spectra of the observed low- and high-mass YSOs. In the case of the high-mass 
sources, the YSOs are disposed as in Fig. lA.2l The line wings of the 3 12 - 3o3 water transition for the high-mass sources have been 
fitted with a Gaussian profile, subtracted and the residuals plotted to isolate the C'^O 7=10-9 emission line. 
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Fig.B.l. Comparison between high-7 HIFI and low-7 JCMT 
spectra for the low-mass source Ser SMMl. '^CO 7=10-9 spec- 
tra (left-top) and C'**0 7=9-8 line (right-top) observed with 
HIFI and for the '^CO 7=3-2 and C'^^O 7=3-2 lines (left-bottom 
and right-bottom respectively) observed with JCMT. The spectra 
have been resampled to 0.27 km s"' and shifted to zero velocity. 
The green line indicates the baseline subtraction. 
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Fig. B.2. Same as Fig. IB.ll but for the intermediate-mass source 
NGC2071. 
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Fig. B.3. Same as Fig. IB.ll but for the high-mass source W3- 
IRS5. The pink Gaussian profile represents the broad velocity 
component identified in the '^CO 7=10-9 line. 
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Table B.l. Values of the FWHM for the '^CO /=3-2 broad ve- 
locity component and the FWHM and integrated intensity for the 
C^^O 7=3-2 spectra. 
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Notes. No broad velocity component identified in the '^CO 7=3-2 
convolved central spectrum. 
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Fig.B.4. '^CO 7=3-2 spectra for low-, intermediate and high-mass YSOs. observed with the JCMT. The green line represents the 
baseline level and the pink Gaussian profile the broad velocity component for those sources for which a 2 Gaussian decomposition 
has been performed. All the spectra have been shifted to zero velocity. The numbers indicate where the spectra have been scaled for 
greater visibility. 



I. San Jose-Garci'a et al.: HIFI CO observations of YSOs: from low- to high-mass. Online Material p 14 



Low— mass: Class 



Low— mass: Class I 



1.2 

- I 
0.8 - 

0.4 - 



IRAS 


2A 




xO.7 - 







B335 









x0,6 



; SMMl 


I . 

xO.5 - 






; SMM3 


xO.5 - 

1 , 



; L723- 

: , r^n^r 


MM 




; L4B3 





L1551 IRS5 



srv|^^-^^iv/-i^ 



Oph IRS63 



' TMRl 

: J 


x2 ' 






' GSS30 


IRSl 




xO.7 


^^^^^ 






L1157 

. r 


x2 





0.8 



L1448 



0.4 - 

Lii^y^ 



IRAS4i 


3 

xO.6 








BHR71 


xO.3 - 





; SMM4 


xO.5 





3^ 



r Elias 


29 : 





^^^^ 



-20 20 



CedllO IRS4 



r L14B9 


x3 - 







-20 20 



-20 20 -20 20 -20 20 

V (km/s) V (km/s) 



-20 20 -20 20 

V (km/s) 



Intermediate— mass 



L1641 S3MMS1 



cq 6 
S 



NGC7129 FIRS2 




-20 



20 -2" ,0 , V 2" -20 20 

V (km/s) 




Hieli- mass 



NGC6334-I(N) W43-MM1 



W33A 







029.96- 


O.OZ 
xO.7 - 





G5.a9-0 


39 

I xO.5 _ 





G34.26+f 


0.15 

1 xO.4 _ 






W51N-el 











NGC753B-IRS1 



-20 20 -20 20 -20 20 -20 20 -20 20 



V (km/s) 



Fig. B.5. Same as Fig. |B3]but for the C 



O J— 3-2 spectra from the low-, intermediate- and high-mass YSOs. 



